Energy separation of single-particle and continuum states in a S = 1/2 weakly-coupled 

chains antiferromagnet. 
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Inelastic neutron scattering is used to study transverse-polarized magnetic excitations in the quasi- 
one-dimensional 5 = 1/2 antiferromagnet BaCu2Si2C>7, where the saturation value for the Neel 
order parameter is mo = 0.12 //b per spin. At low energies the spectrum is totally dominated by 
resolution-limited spin wave-like excitations. An excitation continuum sets in above a well-defined 
threshold frequency. Experimental results are discussed in the context of current theories for weakly- 
interacting quantum half-integer spin chains. 
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One of the outstanding problems in quantum mag- 
netism today is that of dimensional crossover in systems 
composed of weakly-coupled spin chains. In particular, 
for 5=1/2 Heisenberg antiferromagnets (AFs), the in- 
terplay between classical dynamics in Neel-ordered three- 
dimensional materials (single-particle order-parameter 
excitations, or spin waves) and critical dynamics of 
the quantum-disordered 5=1/2 one-dimensional (ID) 
model (2-spinon excitation continuum is not fully 

understood. In quasi-lD 5 = 1/2 systems, for ar- 
bitrary weak inter-chain interactions, long-range order 
is restored. However, according to recent quantum- 
mechanical calculations |3[^], as long as the 3D order 
parameter remains small, the excitation spectrum has a 
unique dual nature. At low frequencies the diffuse con- 
tinuum is "cleaned up" and replaced by long-lived single- 
particle "magnon" excitations. What remains of the crit- 
ical dynamics in isolated chains, is now seen at higher 
energies as a series of multi-magnon states. The lowest- 
energy contribution is from two-magnon processes. At 
the ID AF zone-center the continuum is thus expected 
to set in at twice the characteristic magnon energy, and 
should be separated from the single-particle part of the 
spectrum by a "second gap" . 

To date, this separation of single-particle and contin- 
uum dynamics has not been clearly observed experimen- 
tally in any quasi-lD material. Among other technical 
obstacles, is the limited choice of suitable model com- 
pounds. Most of what is experimentally known of cou- 
pled 5=1/2 chains comes from neutron scattering stud- 
ies of KCUF3 @-p|. Here the ordered moment is rather 
large, making the spin waves very intense at low ener- 
gies, and the continuum contribution difficult (though 
not impossible) to isolate. Materials like Sr2Cu03 ||] 
or SrCu02 Jl^,|ll|, on the other hand, are very good 
ID systems, and it is the spin waves that are hard to 



identify. Below we report inelastic neutron scattering 
studies of spin dynamics in the 5=1/2 quasi-lD an- 
tiferromagnet BaCu2Si2C>7. Using neutron setups with 
complimentary resolution characteristics, we clearly ob- 
serve the frequency-separation of single-particle and con- 
tinuum contributions. In doing so, we pay special at- 
tention to the polarization of magnetic excitations, and 
compare the results to existing theoretical models. 

The recently characterized BaCu2Si2C>7 (orthorhom- 
bic space group Pnma, a = 6.862 A, b = 13.178 A, 
c = 6.897 A) appears to be an ideal model material for 
the present study. Compared to KCUF3 , it has a stronger 
ID character: with an in-chain (along the c-axis) Heisen- 
berg exchange constant J = 279 K, the Neel temperature 
is only Tn = 9.2 K. The structure of the magnetically or- 
dered phase was previously guessed from bulk data and 
an analysis of a single magnetic Bragg reflection |fl^| . As 
part of the present study, we performed a comprehen- 
sive neutron diffraction study of the magnetic structure 
at T = 1.5 K using the D23 single-crystal diffractome- 
ter at Institut Laue-Langevin (ILL) Overall, 35 

inequivalent magnetic Bragg intensities were measured. 
The spin arrangement was found to be exactly as pro- 
posed in Ref. |12|| . The obtained value for the ordered 
moment too = 0-12 /^b is much smaller then in KCUF3 
(too ~ 0.5 /is), yet large enough to expect a measurable 
spin-wave contribution to inelastic scattering. 

Low-energy magnetic excitations in BaCri2Si207 were 
studied in high-resolution inelastic neutron scattering ex- 
periments on a 5 x 5 x 50 mm 3 BaCu2Si207 single-crystal 
sample at the IN14 cold-neutron spectrometer at ILL 
p3| . Measurements were performed in the (0,k,l) scat- 
tering plane with a fixed final neutron energy of 3 meV. 
Typical constant-Q and constant-i? scans collected near 
the ID AF zone-center I = 1 are shown in Fig. 1. Figure 
2 shows a series of constant- Q scans for different mo- 
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mentum transfers perpendicular to the chain-axis. The 
polarization dependence of the magnetic neutron scatter- 
ing cross section ensures that within this range, to a good 
approximation, only fluctuations of spin components per- 
pendicular to the chains (and thus the ordered moment, 
which is parallel to c) are observed. Very similar data 
sets (not shown) were collected using a comparable setup 
on the TASP cold-neutron 3-axis spectrometer at Paul 
Scherrer Institut, Switzerland, for momentum transfers 
along the (h, 0, 1) and (h, h, 1) reciprocal-space rods [ p^[ . 



with J x = -0.463(2) meV, J v = 0.161(1) meV and 
2J 3 = 0.145(1) meV. Simulations based on these values 
are shown in solid lines in Figs. |l|,||. The resulting dis- 
persion relation along the (0, k, 1) direction is shown in 
a solid line in the (Tiuj)-k plane in Fig. |[ Additional 
analysis also gives an upper estimate for the intrinsic 
excitation width: 0.03 meV. These results conclusively 
demonstrate that up to about 4 meV energy transfer the 
magnetic excitation spectrum is fully accounted for by 
the single-particle picture, in the entire range of momen- 
tum transfers perpendicular to the chain-axis. 
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FIG. 1. Typical constant-Q (a,b) and constant-i5 (c) scans 
measured in BaCu2Si207 at T = 1.5 K near the ID AF 
zone-center I — 1 (symbols). Lines are as described in the 
text. 

The main feature seen in all constant-Q scans is a sharp 
asymmetric peak that we attribute to single-magnon ex- 
citations. The intensity onset on the low-energy side is 
very steep, and suggests the peaks are resolution-limited. 
The extended high-energy "tail" is due to a steep disper- 
sion along the chain axis and a non-zero wave vector res- 
olution in that direction. The data were analyzed using a 
single-mode approximation (SMA) cross section, derived 
from the quantum chain-Mean Field (chain-MF) model 
j|. The dispersion measured along (0, k, 1), (h, 0, 1) and 
(h, h, 1), we found that at least three inter-chain ex- 
change constants are needed, between nearest-neighbor 
spins along the (1,0,0), (0,1,0) and (1,1,0) real-space 
directions, that we denote as J Xl J y and J3, respectively. 
The SMA cross section was convoluted with the calcu- 
lated experimental resolution function and used to fit all 
measured scans to determine J x , J y and J3. The in- 
chain coupling constant J was fixed at J — 24.1 meV, 
as previously deduced from bulk susceptibility curves 
p2[. An excellent global fit to all scans is obtained 
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FIG. 2. A series of constant-Q scans measured in 
BaCu2Si207 at T — 1.5 K for different momentum trans- 
fers perpendicular to the chain axis. Lines are as described in 
the text. 

Extending the cold-neutron study to higher energy 
transfers was hindered by certain geometrical constrains 
imposed by the design of IN14 spectrometer. Instead, the 
intermediate-energy spectrum was studied using the IN22 
thermal-neutron spectrometer at ILL using 14.7 meV 
final-energy neutrons. Most of the data were collected 
in constant- E scans along the chain-axis, in the vicinity 
of the (0, 0, 3) ID AF zone-center. The background orig- 
inating from the empty Al sample-holder was measured 
separately for each scan. Typical background-subtracted 
data sets are presented in Fig. ||. None of these show 
the distinct 2-peak structure as at 3 meV energy transfer 
(Fig. lc). To verify that this is not a resolution effect, 
we performed a least-squares fit of the SMA cross sec- 
tion, convoluted with the calculated resolution, to each 
measured profile. Very poor fits are obtained, as shown 
in dashed lines in Fig. |[ It is clear that if the observed 
intensity was due to a single mode, two separate peaks 
would have been easily resolved. For lack of a conve- 
nient analytical result for interacting chains, we fitted 
the measured scans to the Miiller-ansatz continuum form 
|l4j], known to work well for isolated chains. The agree- 
ment in this case is much better (solid lines in Fig. |^). 
For contrast, the dashed line in Fig. |l|c shows the best 
Miiller-ansatz fit to the 3 meV high-resolution const--E 
scan. From this comparative analysis we conclude that, 
unlike at low energies, a large contribution to the ob- 
served scattering at high energies must come from con- 
tinuum excitations. 

Additional evidence was obtained in analyzing the 
measured intensities. In the chain-MF model the in- 
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tensity of magnon scattering has the same l/u> depen- 
dence as for classical antiferromagnetic spin waves, in 
full agreement with the bulk of our cold-neutron data. 
Having confirmed the single-mode picture at low ener- 
gies, we exploited this scaling relation to estimate spin 
wave intensities at higher energy transfers. A normaliza- 
tion factor was obtained by fitting the SMA cross sec- 
tion to the intensity measured using the thermal setup 
at q = (0, 0, 3) between 2 and 3 meV transfer, where, ac- 
cording to our cold- neutron results, the SMA picture is 
still valid. The thus calculated single-mode contributions 
are represented by shaded areas in Fig. ||. 
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FIG. 3. Constant-i? scans measured in BaCu2Si207 at 
T = 1.5 K using a thermal-neutron setup. Lines are fits 
to the data, as discussed in the text. Shaded areas are the 
extrapolated single-mode contributions. 

Clearly, magnon excitations account for only a minor 
fraction of the dynamic susceptibility at high frequencies. 
To quantify this behavior, in Fig. [|a we plot the mea- 
sured g z -integrated intensity I z , scaled by energy trans- 
fer Tluj, as a function of the latter. The advantage of 
this integration procedure that the actual wave-vector 
dependence of continuum scattering becomes unimpor- 
tant. The expected spin wave contribution is practically 
constant above 3 meV in this plot (Fig. |]a, dashed line). 
This follows from that in a fixed-final energy configura- 
tion one directly measures S(q, u>) without any additional 
energy-dependent scaling factors (A/2 effects in the mon- 
itor are expected to be less than 5%, thanks to the use 



of a thermal guide). In contrast, experimentally, I z u> in- 
creases substantially with frequency. 

The continuum fraction can be isolated by subtracting 
the extrapolated SMA part from the experimental data. 
The result is shown in Fig. ^b. The energy resolution 
being about 1.5 meV, we conclude that the continuum 
sets in between 3.5 and 5.5 meV energy transfer. We 
have verified that this result is quite robust, and remains 
valid even if we assume that there is a 50% systematic er- 
ror in the estimate of the single- mode contribution. From 
Fig. |]b we see that above the threshold the (^-integrated 
continuum intensity is only slightly energy-independent. 
Looking back at the cold-neutron measurements, we find 
that they actually do contain direct evidence of contin- 
uum scattering as well. Indeed, the increase of intensity 
seen around 4.5 meV in Fig. [l] is not accounted for by 
the SMA picture, and thus represents the onset of the 
continuum, in full agreement with the thermal-neutron 
results. 
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FIG. 4. (a) Measured g z -integrated intensity scaled by en- 
ergy transfer (symbols). The dashed line is the extrapolated 
single-mode contribution, (b) Estimated ^-integrated inten- 
sity of continuum excitations in BaCu2Si2 07 (symbols). The 
dash-dot line is the calculated spin wave theoretical 3-magnon 
spectrum. In both panels the solid line is a guide for the eye. 

To understand the observed long-lived and continuum 
excitations, we shall recall the simple yet profound phys- 
ical picture provided by the quantum chain-MF theory 
HQ. For an isolated chain in a staggered field, the 
magnons have an energy gap A that scales as H^ 3 
p5|,p|p6|. In the coupled-chains case, magnon dispersion 
perpendicular to the chain axis leads to a softening of the 
gap at the 3D zone-center for transverse-polarized exci- 
tations. These gapless modes correspond to conventional 
spin waves in the classical spin wave theory (SWT), and 
are the long-lived excitations that we see in BaCu2Si207 
at low energies. Excitations at exactly A can be ob- 
served in points of reciprocal space where inter-chain- 
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interactions cancel out at the MF level. In BaCu2Si20y 
this occurs at q = (0.5,0.5,1). Using the measured ex- 
change constants we can estimate A = 2.37 meV. It is 
straightforward to verify that the experimental values for 
A, J, mo and Tn satisfy the relations predicted by the 
chain-MF model || remarkably well. 

A unique quantum-mechanical feature of a half-integer 
Heisenberg AF chain in a staggered field, totally ab- 
sent in conventional SWT, is a longitudinal mode, i.e., 
a magnon excitation polarized parallel to the direction 
of static staggered magnetization. The existence of such 
an excitation was recently confirmed experimentally in 
KC11F3 0. Focusing on transverse spin correlations in 
the present study, we do not directly observe this branch 
in BaCu2Si207. It does, nonetheless, play an important 
role in the observed transverse continuum. In the chain- 
MF and field-theoretical models for weakly-interacting 
chains, low-energy continuum excitations are seen as 
multi- magnon states B. The existence of a longitudi- 
nal mode allows iwo-particle excitations in the transverse 
channel, composed of one longitudinal and one transverse 
magnon. The threshold for such excitations at the ID AF 
zone-center is at twice the characteristic magnon energy: 
A c = 2A. This prediction is consistent with our experi- 
mental results for BaCu2Si207, where A c is between 3.5 
and 5.5 meV, and 2A = 4.8 meV. 

An alternative model that predicts multi-magnon con- 
tinua is SWT with kinematic corrections, which explic- 
itly takes into account non-linearities in the spin wave 
Hamiltonian. Continuum excitations observed in KCUF3 
are remarkably well described by this model In 
BaCu2Si20y, however, the ordered moment is much 
smaller, and the ground state resembles the Neel phase, 
i.e., the starting point of SWT calculations, consider- 
ably less. For this material SWT completely fails at the 
quantitative level. For example, the spin wave correction 
to sublattice magnetization in BaCu2Si20y that we cal- 
culated using the known dispersion bandwidths is well 
over 100%. For the transverse-polarized excitation con- 
tinuum, the inadequacy of SWT becomes even more ap- 
parent. Indeed, the longitudinal mode being absent in 
this model, two-magnon transverse excitations are pro- 
hibited by selection rules. The lowest-energy contribu- 
tion is from three-magnon states. The lower bound of the 
continuum should thus be at roughly 3 A. The 3-magnon 
SWT contribution can be exactly evaluated, provided 
the spin wave dispersion is known. Following the recipe 
given in Ref. , the spin wave bandwidths measured in 
BaCu2Si207 were used to calculate the 3-magnon profile 
at q = (0,0,3). The result (arbitrary scaling) is shown 
in Fig. |^b in a dash-dot line. The calculated continuum 
onset is clearly at a larger energy, by 2 to 3 meV, than 
observed experimentally. 

In summary, we find that the transverse excitation 
spectrum in BaCu2Si2C>7 is divided into two well-defined 
regions. Below 4 meV the weight is consolidated into 



long-lived single-particle spin wave-like excitations, and 
no sign of multi-particle continuum is found. These ex- 
citations carry only a small fraction of the total spectral 
weight at higher energies, where an excitation continuum 
becomes the dominant contribution. While further stud- 
ies of the lower continuum bound are needed, the present 
results suggest a threshold between 3.5 and 5.5 meV. For 
this strongly ID material the conventional spin wave- 
theoretical description becomes inadequate, even if kine- 
matic corrections are taken into account. In contrast, the 
chain-MF model, based on quantum-mechanical proper- 
ties of individual chains, is in very good agreement with 
experiment. 
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